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The dynamics of the one-dimensional g-state Potts model, in the zero temperature limit, can be 
formulated through the motion of random walkers which either annihilate (A + A — » 0) or coalesce 
(A + A — > A) with a g-dependent probability. We consider all of the walkers in this model to be 
mutually infectious. Whenever two walkers meet, they experience mutual contamination. Walkers 
which avoid an encounter with another random walker up to time t remain uninfected. The fraction 
of uninfected walkers is known to obey a power law decay U(t) ~ t~^ q \ with a nontrivial exponent 
<t>{q) [C. Monthus, Phys. Rev. E 54, 4844 (1996); S. N. Majumdar and S. J. Cornell, ibid. 57, 
3757 (1998)]. We probe the numerical values of (j>(q) to a higher degree of accuracy than previous 
simulations and relate the exponent <f>(q) to the persistence exponent 8(q) [B. Derrida, V. Hakim 
and V. Pasquier, Phys. Rev. Lett. 75, 751 (1995)], through the r elation (j>(q) = , y(q)6 (q) where 7 is 
a new exponent introduced in [S. J. O'Donoghue and A. J. Bray, cond-mat/01 11133] ]. Our study is 
extended to include the coupled diffusion-limited reaction A + A— * B, B + B — > Ain one-dimension 
with equal initial densities of A and B particles. We find that the density of walkers decays in this 



model as p(t) ~ t 



-1/2 



The fraction of sites unvisited by either an A or a B particle is found to 



obey a power law, P(t) ~ t with 9 ~ 1.33. We discuss these exponents within the context of 
the q-state Potts model and present numerical evidence that the fraction of walkers which remain 
uninfected decays as U(t) ~ t~^, where <f> ~ 1.13 when infection occurs between like particles only, 
and (j> ~ 1.93 when we also include cross-species contamination. We find that the relation between 
<f> and 8 in this model can also be characterized by an exponent 7, where similarly, (j> — yO, 



X 



I. INTRODUCTION 



When a system is quenched from a homogeneous high- 
temperature disordered state into a low-temperature 
regime, well defined domains of equilibrium ordered 
phases form randomly and then grow with time 
("coarsen") in a self-similar way, until the domain size 
becomes comparable to the size of the system Q. A 
dynamic scaling hypothesis suggests that, at late times, 
the evolution of the system is characterized by a sin- 
gle length scale L(t) that represents a typical linear size 
of the domains. It is well established, at least for sys- 
tems with a scalar order parameter, that L(t) ~ t n with 
n = 1/2 for non-conserved dynamics and n — 1/3 for 
conserved dynamics Q. The Ising spin model evolving 
with Glauber spin flip dynamics pj demonstrates behav- 
ior characteristic of the former while evolution according 
to Kawasaki spin dynamics ]3| exemplifies the latter. A 
generalization of the Glauber-Ising model is the g-state 
Potts model which has q distinct, but equivalent, 

ordered phases. Experimental realizations arc known for 
q = 2 i.e. the Ising model and for q = 3,4, 00 |4[. The 
q = 00 case describes several cellular structures [S, e.g . 
polycrystal s M] , soap froths Jll[], magnetic bubbles fj^j 
and foams Jl3|. 

Such coarsening systems are amongst those which have 
received considerable attention in recent years with re- 
gard to their persistence properties [|l3-48|. The persis- 
tence probability P(t) is simply the probability that a 
given stochastic variable (f>[x,t), with zero mean, does 
not change sign in the time interval [0, t]. Theoretical 



and computational studies of persistence include spin 
systems in one p^|,|T6| and higher JT7|-pl| dimensions, 
diffusion fields |22|,|23|, fluctuating interfaces pi] , phase- 
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ordering dynamics |p5|-[27[ , Lotka-Volterra models 
and reaction-diffusion systems [p0[-p8[. Experimental 
studies include the coarsening dynamics of breath fig- 
ures [|9| , foams Jl3| , soap froths @,|dj , twisted nematic 
liquid crystals |42| ], and one-dimensional gas diffusion 
|43| . Persistence in non-equilibrium critical phenomena 
has also been studied in the context of the global order 
parameter M(t) (e.g. the total magnetization of a fer- 
romagnet), regarded as a stochastic process p^ , p5[ . In 
many systems of physical interest, the persistence decays 
algebraically according to P(t) ~ t~ 9 , where 9 is the per- 
sistence exponent and is, in general, nontrivial. The non- 
triviality of 6 emerges as a consequence of the coupling of 
the field 4>{x, t) to its neighbors, since such coupling im- 
plies that the stochastic process at a fixed point in space 
and time is non-Markovian. For the g-state Potts model 
at T = 0, the fraction of spins which have never flipped 
up to time t (i.e that fraction that have persisted in their 
original phase), has been observed [}l5|Jl7| -|l9[| to obey a 
power law decay P(t) ~ t~ et - q \ where 9(q) has been ob- 
tained exactly, in one-dimension (ID), by Derrida et al. 
p6| for arbitrary q, 
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There is a very direct way |49|-|5l|] of relating the Glauber 
dynamics of the ID g-state Potts model at zero temper- 
ature to reaction-diffusion models [B2l,B3|. Consider un- 
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correlated initial conditions where each of the q phases is 
present with equal density l/q. During any time interval 
dt, each spin Si(t) has a probability dt/2 of becoming 
equal to its right neighbor, dt/2 of becoming equal to 
its left neighbor and a probability 1 — dt of retaining its 
own value. The domain walls therefore perform random 
walks. Upon contact, two domain walls react instanta- 
neously, either annihilating A + A — > with probabil- 
ity l/(q — 1) or coalescing A + A — > A with probability 
(q— 2)/(q — 1), these numbers being the probabilities that 
the states on the far sides of the walkers are the same 
(annihilation) or different (coalescence). This reaction- 
diffusion model, together with the fact that in the initial 
condition each bond is occupied by a domain wall with 
probability (q— l)/q, is completely equivalent to the spin 
problem. The fraction of spins which have never flipped 
up to time t is then simply the fraction of sites which have 
never been visited by a random walker. The probabilistic 
algorithm for implementing the Potts model through the 
annihilation or coalescence of random walkers allows q to 
be a real number but is restricted to q > 2. However, 
an equivalent Isingspin representation of the Potts per- 
sistence problem M] permits a study |37| of Q(q) in the 
regime 1 < q < 2. 

In the present work we study the fraction of random 
walkers which have never encountered another walker, 
using the reaction-diffusion representation of the ID 
Glauber-Potts model at zero temperature. To facilitate 
our discussion, we consider all the random walkers at 
time t = to be uninfected. All the walkers are, how- 
ever, considered to be infectious, so that for t > any 
contact between walkers (assuming they survive the en- 
counter) leads to mutual contamination. Walkers that 
avoid all contact remain uninfected. Our goal, therefore, 
is to address the fraction of uninfected walkers U(t) up 
to time t. 

This problem has been addressed in some detail by a 
number of different authors |2^,^,^], although under 
the guise of a different interpretation. The general ap- 
proach used involves monitoring the motion of a tagged 
test particle, released in the system at time t = 0. If the 
tagged particle diffuses with diffusion constant D' and 
the other particles in the system diffuse with diffusion 
constant D, the tagged particle is viewed as an external 
impurity diffusing through a homogeneous background. 
One is then interested in addressing the survival prob- 
ability of the tagged particle i.e., the probability is has 
not been absorbed by another particle in the system. Of 
course in our terminology, this corresponds to the prob- 
ability that the tagged particle remains uninfected. It is 
well established |34|,[47J that U(t) decays according to a 
power law U(t) ~ t^^ 9,c \ where c = D'/D is the ra- 
tio of the two diffusion constants in the problem. When 
D' = 0, corresponding to a static impurity, the model 
clearly maps onto the problem of persistent spins and in 
this case, <f>(q, 0) = 9{q). The exponent <f>(q, c) is therefore 
considered a generalization of the standard persistence 
exponent 6(q), characterizing the survival probability of 



a mobile particle under the coarsening dynamics of the 
Potts model. A recent study of the survival probability of 
a mobile particle, moving according to either determinis- 
tic or stochastic rules, in given in |47|] , where the external 
fluctuating field takes the form of either the solution to 
the diffusion equation, the coarsening dynamics of the 
q-state Potts model or spatially uncorrelated Brownian 
signals. 

In the limit q — > oo the Potts model is equivalent to the 
A + A — > A single species coalescence process. The dy- 
namics of this system are particularly simple, since each 
particle sees only the "cage" formed by its two nearest 
neighbors. It is immaterial that these neighbors may 
coalesce with more distant particles. Consequently, the 
kinetics of each particle involves only itself and its two 
nearest neighbors and one is therefore able to calculate 
0(oo,c) exactly 



(oo,c) 



2cos~ 



l+c 



(2) 



(Note that for c = 0, this of course reduces to the static 
persistence exponent, 0(oo, 0) = 6(oo) = 1.) In the limit 
of the Ising model (q — 2), there is, however, no cor- 
responding exact result for general c. An exact result 
is only known for c = i.e., 0(2,0) = 9(2) = 3/8. A 
mean-field Smoluchowski ]2(| approach predicts for q = 2 
and general c that, 0(2, c) — + c)/8, in good agree- 
ment with numerical simulations |2^| . However, when ex- 
tended to the large q limit, this Smoluchowski approach 
differs substantially from the exact q — > oo result Eq. (||). 
More recently, however, a perturbation theory developed 
by Monthus Q] has evaluated the exponent <fi(q, c) to 
first order in q — 1 for arbitrary c and at first order in c 
for arbitrary q. In this paper, we will not study the expo- 
nent <f)(q, c) in complete generality since we are primarily 
interested in the special case where all of the particles 
are equally mobile (D 1 = D) and henceforth we therefore 
write 4>(q, 1) = 4>{q). 

Within the context of the current work, where we are 
interested in the case c = 1, the random walkers repre- 
sent the motion of domain walls undergoing zero temper- 
ature coarsening. The fraction of uncollided domain walls 
(i.e., uninfected walkers) has previously been interpreted 
as "domain wall persistence" or equally, the probability 
that two adjacent domains survive. Under this interpre- 
tation, Majumdar et al. have studied the quantity 
U(t) numerically for various q. They observed a power 
law decay, U(t) ~ t^^^ in agreement with the results 
of Monthus 134] and obtained a spectrum of values for 
4>(q): 1/2 < 4>(q) < 3/2, where 2 < q < oo. In this pa- 
per, we extend and improve these numerical results by 
presenting (Sec. II) more accurate data than previously 
achieved EtJ and exploring a larger range of q. We also 
relate the value of <fi(q) to the static persistent exponent 
0(q) via <p(q) — j{q)0(q), where 7 is a new exponent that 
we recently introduced in [p7[. 
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We additionally consider (Sec. Ill) the coupled 
diffusion- limited reaction A + A-^ B, B + B — > Am one- 
dimension with equal initial densities, pA (0) = ps(0). 
We find that the density of particles decays in this model 
according to p ~ i -1 / 2 , independent of the initial walker 
density, characteristic of "Potts like" behavior [pi]] . The 
fraction of sites unvisited by either an A or a B par- 
ticle decays as Pit) ~ t~ e with ~ 1.33. We discuss 
these exponents within the context of the g-state Potts 
model and observe that the fraction of walkers which re- 
main uninfected by their own species decays according 
to U(t) ~ t~^, where cf) ~ 1.13. We also study the frac- 
tion of walkers which remain uninfected by both species 
and find that this too obey a power law, with <fi — 1.93. 
In analogy to Sec. II we express the exponent <fi is terms 
of the static persistence exponent 8 and similarly obtain 
4> = 7$. We conclude in Sec. IV with a summary of the 
results and a discussion of some open questions. 



II. THE POTTS MODEL 



We investigate numerically the fraction of uninfected 
walkers U(i) in the ID g-state Potts model, at zero tem- 
perature, as a function of q. Our simulations are per- 
formed on a ID lattice of size 10 7 with periodic boundary 
conditions. At t = a random walker is placed at each 
lattice site, so that p(0) = 1. We choose this initial high 
density of walkers to accelerate the system's evolution 
into the asymptotic regime. Our model is updated using 
the direct method |55| i.e. at each computational step, a 
particle is picked at random and moved with probability 
D = 1/2 to a neighboring site, where D is the diffu- 
sion constant. If the destination site is occupied, the two 
particles either annihilate, with probability l/(g — 1), or 
coalesce, to become a single particle, with probability 
iq — 2)/iq— 1), in accordance with the reaction-diffusion 
dynamics of the Potts model detailed in Sec. I. For each 
move of a particle, time t is increased by dt — 1/N, where 
N is the current number of particles in the system. Such 
sequential dynamics can be chosen without loss of gen- 
erality as parallel dynamics exhibit similar asymptotic 
behavior [p6| . Our simulations are performed for 5000 
time steps and our results are averaged over 100 runs. 

Fig. |l| clearly shows that the fraction of uninfected 
walkers decays according to a power-law, U(t) ~ t~^ q \ 
In Table I we present our numerical values for <f>iq) in 
the range 2 < q < oo. These values were obtained by 
performing a linear regression on log-log curves such as 
those presented in Fig. |l|. The regression was taken in 
the range 10 < t < 1000. This range was chosen to 
avoid initial transients and to avoid statistical fluctua- 
tions between different runs which become prominent as 
t increases, especially for large q. In Fig. |2] we plot <j> as a 
function of q, where we choose to represent the data on 
a linear-log plot merely for clarity of presentation. 
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FIG. 1. Log-log plot of the fraction of uninfected walkers 
in the ID, T = 0, Glauber-Potts model as a function of time, 
for various q. 



q 





2 


0.5006(7) 


4 


0.867(5) 


8 


1.121(5) 


16 


1.283(6) 


25 


1.352(6) 


32 


1.380(7) 


50 


1.419(7) 


64 


1.435(8) 


100 


1.461(8) 


128 


1.464(8) 


256 


1.481(9) 


oo 


1.495(9) 



TABLE I. Numerical values of the exponent <j>{q) in the ID, 
T = 0, Glauber- Potts model where U(t) ~ t'^, for various 
values of q. 

The case q = 2, which corresponds to the Ising model, 
reduces to the A + A — > reaction-diffusion process, 
where the density of uninfected walkers is clearly equal 
to the particle density. It is well known |54| that in the 
g-state Potts model the particle density obeys a power- 
law decay, pit) ~ i -1 / 2 , which is consistent with our 
value of 0(2) = 0.5006(7). For the case q = oo, one 
can obtain the exact value of <p(oo) from Eq. (0) i.e. 
</>(oo, 1) = 3/2. This is again consistent with our nu- 
merical result, <fi(oo) = 1.495(9). These limiting values, 
along with our numerics, suggest that the fraction of un- 
infected walkers decays according to U(t) ~ t^^^ where 
1/2 < 4>iq) < 3/2, for 2 < q < oo, consistent with the 
results of Monthus Q and Majumdar et al. |47). Al- 
though the values of (j>iq) = 1/2 and (/>(oo) = 3/2 at the 
limit of the q spectrum are understood, it was pointed 
out by Monthus that the formalism that led to the 
exact determination of the static persistence exponent 
9(q) |Q cannot be easily applied to compute the exact 
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value of 4>(q,c) for general c, including the limit c — I. 
The exact calculation of 4>(q) therefore remains an open 
problem. 




FIG. 2. Linear-log plot of as a function of q for the ID, 
T = 0, Glauber-Potts model. 

Recently, we have explored the connection between 
uninfected walkers and unvisited sites in one-dimension 
within the context of a system of non-interacting ran- 
domly diffusing particles A% < — ► 0A, and the A + B -> 
diffusion- limited reaction J57J. We found that in both of 
these models, 

U(t) ~ [P(£)] 7 , with 7 ~ 1.39. (3) 

While this relationship is obeyed reasonably accurately 
over the range of times accessible to our simulations, 
there is evidence, for both processes, that 7 approaches a 
smaller value at late times |57| . Clearly, we can write an 
analogous relationship for the g-state Potts model, where 
U(t) ~ t-^\ P(t) ~ t- 6 ^ and 7 is now some function 
of q, i.e. t-'t'W = t-f^ e( i\ so that 

<P{q) = l{q)0{q) ■ (4) 

We plot <p(q) against 9(q) in Fig. 0, where our values for 
4>(q) have been taken from Table I and the corresponding 
values of 0{q) have been calculated exactly from Eq. ([I]). 
The point at the origin corresponds to q = 1, where </> 
and 6 both vanish. Fig. || indicates that 0(g) increases 
monotonically as a function of 6{q). Therefore, given that 
we know the exact values of 9{q) and <p(q) for q = 2, 00, 
and given the convex nature of the plot in Fig. |3j, we can 
use Eq. (Q) to determine precise bounds on 7(g) in this 
range: 

4/3 < 7(9) < 3/2, ge[2,oo]. (5) 

Thus while 4>{q) and 6(q) are separately quite strongly 
dependent on g, their ratio 7(g) depends only weakly on 

q- 
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FIG. 3. The exponent <f> plotted against the exponent 8, 
for various q for the ID, T=0, Glauber-Potts Model. 



III. THE EXCHANGE MODEL 

The A + A^B, B + B^A diffusion- limited re- 
action offers a curious combination of birth and death 
processes. For brevity, we dub this model the "exchange 
model" . The particles in the exchange model diffuse ran- 
domly. Whenever a particle encounters a walker belong- 
ing to its own species, the two particles annihilate in- 
stantaneously, to create a particle belonging to the other 
species at the point of annihilation. For the purposes 
of this study, we allow lattice sites to contain multi- 
ple particles, enabling two particles of differing species 
to cohabit a single lattice site. This feature of our 
model ensures that, in one-dimension, the lattice struc- 
ture ...ABABABABABAB... continues to evolve. Con- 
fining our system to a maximum of one particle per site 
would otherwise mean that the preceding configuration 
would be stable for all time t. We focus our attention on 
the ID exchange model with equal initial densities of A 
and B particles. 

Our numerical simulations are performed on a one- 
dimensional lattice with periodic boundary conditions. 
At t = equal numbers of A and B particles are ran- 
domly distributed on the lattice with a maximum of one 
particle per site. Our system then evolves according to 
the dynamics described above. The system is updated in 
the same manner as for the g-state Potts model described 
in Sec. II. 



A. Particle Density 

The first question we address is the particle decay. The 
unimolecular reactions defining the exchange model are 
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characteristic of the g-state Potts model and therefore we 
can expect the particle density to exhibit similar asymp- 
totic behavior. Indeed, for every pair of reactions be- 
tween two A and two B particles, a single A and a single 
B are produced, which is the same result as for two q = oo 
Potts models, A + A — ► A, B + B — > B operating on a 
single lattice. We therefore expect the particle density to 
decay in the exchange model according to pAif) ~ Ct~ a , 
with a = 1/2, independent of the initial particle density 

f4j. This is indeed confirmed by our simulations. In Fig. 
we present our results for the particle density for three 
different initial values of Pa{0) (recall ps(0) = pa(0)). 
The simulations were performed on a lattice of size 10 6 
for 10 5 time steps and the results averaged over 100 runs. 
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FIG. 4. Log-log plot of the particle density as a function of 
time in the ID exchange model. 

In Table II we present our numerical values for the par- 
ticle density decay exponent a and amplitude C, where 
PA(t) ~ Ct~ a . These values were obtained by perform- 
ing a linear regression on the the plots in Fig. ^ in the 
asymptotic regime, 10 4 < t < 10 5 . 



Pa(0) 


a 


C 


0.5 
0.25 
0.1 


0.5007(9) 
0.5007(9) 
0.5000(5) 


0.419(5) 
0.419(5) 
0.416(9) 



TABLE II. Numerical values for the particle density decay 
exponent a and the amplitude C in the ID exchange model, 
where pA(t) ~ Ct~ a , for different initial densities Pa(0). 

For the ID g-state Potts model, the amplitude C is 
g-dependent p%, 



C{q) 



1 1 



q V2ttD 



(6) 



Is there perhaps a value of q which corresponds to the 
exchange model? A value of q = 4 in Eq. (||), with 



D = 1/2 returns a value C(4) = 0.423, in good agree- 
ment with the values of C recorded in Table II. Notice 
that this value lies exactly midway between the min- 
imum value, C(2) = 0.282 and the maximum value, 
C(oo) = 0.564, set by the g-state Potts model. Our nu- 
merical results therefore suggest, that in terms of the 
particle density, the two reactions which define the ex- 
change model (A + A ~ > B,B + B ^ A) are equivalent, 
in some loose sense, to two independent g = 4 Potts mod- 
els, though we would not wish to stretch this comparison 
too far. 



B. Fraction of Unvisited Sites 



We now turn our attention to the persistence proper- 
ties of the exchange model. Here, we are interested in 
the fraction of sites, P(t), that have never been visited 
by either an A or a B particle up to time t. An obvious 
starting point is to see if P(t) obeys a power-law decay 
with an exponent equal to that returned in the case of 
two q = 4 Potts Models i.e. P(t) ~ i- 2x9 (<7= 4 ) where, us- 
ing Eq. (]]]), 2 x 9{q = 4) = 1.263. Our numerical simula- 
tions do indeed indicate a power-law decay, P(t) ~ t~ e " . 
However, the decay exponent 9 ex ~ 1.33, in disagreement 
with what one might have naively expected from a study 
of the particle decay. 

Our simulations are performed on a lattice of size 10 6 
for 5000 time steps and the results averaged over 200 
runs. In Fig. [s] we present a log-log plot of P(t). In Table 
III we show the results of a linear regression performed 
on the curves in Fig. 0, to obtain values of 8 ex . The re- 
gression was performed in the regime 500 < t < 5000, to 
avoid initial transients. 
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FIG. 5. Log-log plot of the fraction of unvisited sites as a 
function of time in the ID exchange model. 



5 



Pa(0) 




0.2 
0.3 


1.330(7) 
1.333(5) 



TABLE III. Numerical values for the persistence decay ex- 
ponent 9 ex in the ID exchange model, where P(t) ~ t -0 " , 
for initial densities Pa(0). 



C. Fraction of Uninfected Walkers 



Finally, we consider the fraction of uninfected walkers 
in the exchange model. We focus, without loss of gen- 
erality, on just the fraction of uninfected A particles. In 
this model, we are only interested in original A particles. 
Those A particles which have been created through the 
reaction B + B — > A after t = 0, are also therefore con- 
sidered to be infected. There are two distinct cases to 
consider: (i) the fraction of A particles which remain un- 
infected by only their own species and (ii) the fraction of 
A particles which have avoided infection by both types of 
species. Our approach is wholly numerical. We address 
the former case first. 

Fig. ^| indicates that the fraction of uninfected A walk- 
ers obeys a power-law decay according to U(t) ~ t~^ A , 
where the subscript A on the exponent denotes infection 
by only A particles (case (i)). The results for pa(0) = 1 
were generated from 500 runs on a lattice of size 5 x 10 5 
over 5000 time steps. The results for pa(0) = 0.5 were 
generated using a lattice of size 10 6 , run for 10 4 time 
steps and also averaged over 500 runs. In Table IV, we 
present our numerical values of 4>a for two different ini- 
tial starting densities. The values of the exponent were 
obtained by performing a linear regression on the curves 
in Fig. | in the regime 100 < t < 5000. 
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FIG. 6. Log-log plot of the fraction of A walkers uninfected 
by A particles in the ID exchange model as a function of time. 



MO) 


4>A 


1.0 
0.5 


1.133(5) 
1.138(8) 



TABLE IV. Numerical values for the exponent <f>A in the 
ID exchange model, where U(t) ~ t~^ A , for initial densities 
Pa(0). 

We now address case (ii), the fraction of A particles 
which have avoided infection by both A and B walkers. 
Our simulations are performed on a lattice of size 10 6 for 
2000 time steps and averaged over 500 runs. In line with 
our previous data, we present our results for U(t) on a 
log-log plot in Fig. ^[ 
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FIG. 7. Log-log plot of the fraction of A walkers uninfected 
by both A and B particles in the ID exchange model as a 
function of time. 

Unfortunately, we do not obtain straight line graphs. 
However, an inspection of individual runs of the data 
(Fig. ||) shows that, in this model, statistical fluctuations 
between different runs become prominent very quickly, 
whereas earlier times are governed by transient effects. 
The dominance of these two regimes renders any deter- 
mination of an exponent particularly difficult. Using a 
lower initial density, with the intention of delaying the 
onset of noisy data only extends the transient regime, 
whereas using a higher initial density to minimize tran- 
sient effects accelerates the onset of large statistical fluc- 
tuations. In view of the power-law behavior we identified 
in case (i) and the general trend of our data, we suggest 
that, with sufficiently good statistics, the fraction of A 
walkers uninfected by either A or B particles would also 
exhibit a power-law decay U(t) ~ t~^ AB , where the sub- 
script on the exponent denotes infection is permitted by 
both species. There is, however, a cleaner method of ex- 
tracting the exponent 4>ab ■ We suggest, in analogy with 
the other models we have studied, that U(t) ~ [P(i)] 7ex 
so that, 

4>AB = lexOex- (7) 
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FIG. 8. Log-log plot of the fraction of A walkers uninfected 
by both A and B particles in the ID exchange model as a 
function of time, for five individual simulation runs. 

If we can establish the validity of the relation U(t) ~ 
[P(t)] 7ex , and determine the value of j ex , we can then 
read off the exponent 4>ab horn Eq. (Q). Our initial 
condition specifies a maximum of one particle per site, 
so J7(0) = 1. While, in the continuum description of 
the model, P(0) = 1 also, on a lattice this becomes 
P(0) = 1 — po- Therefore, to place these two quantities 
on an equal footing in our numerical study, we define, 
p{t) = P(t)/P(0), so that both U(0) and p(0) are equal 
to unity. Note, of course, that in the continuum limit, 
po — > 0, we recover p(t) — > P(t). We present in Fig. |9| a 
log-log plot of U (t) vs. p(t) . 
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FIG. 9. Log-log plot of the fraction of A walkers uninfected 
by both A and B particles, against the fraction of unvisited 
sites, in the ID exchange model. 

The plot in Fig. ^ is quite linear, but does not go through 
the origin, suggesting that in fact U (t) ~ Kp(t) lcx with 



K 1. We note that in the numerical analysis of the 
relationship between noninfectedness and persistence in 
the A% <-> %A and A + B — ► models Q the equivalent 
prefactor K was close enough to unity (K ~ 0.99) to 
tempt the speculation that it might become unity in the 
continuum (low density) limit, though we stress that in 
the latter models there was good evidence that some of 
the simulations had not reached the asymptotic regime. 
In order to evaluate j ex and K for the present model, we 
perform a linear regression on the data in Fig. |^ in the 
region indicated by the arrows, thereby avoiding initial 
transients associated with lattice effects, and the onset of 
statistical fluctuations between different runs of the data. 
We summarize our results for and K in Table V. 



Pa(0) 


lex 


K 


0.2 
0.3 


1.453(7) 
1.456(7) 


0.863(9) 
0.754(9) 



TABLE V. Numerical values of y ex and K in the ID ex- 
change model, where U(t) ~ if[p(t)] 7<!a= , for two initial densi- 
ties Pa(0). 

Note that the value of K for pa(0) — 0.2 is larger than 
that for pa{0) — 0.3, suggesting the possibility that 
K — ► 1 as pa(0) — ► 0, consistent with our results for 
other models J57j , though we have no concrete arguments 
to support this idea. To test our results, we attempt 
to collapse the data for the uninfected fraction and the 
persistence on to a single curve by plotting In U (t) and 
7 ex lnp(t) + ln.fr as functions of time t, in Fig. nG. 
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FIG. 10. In U (t) and "/ex lnp(t) +ln K plotted as a function 
of time, to collapse the data for the uninfected walkers and 
unvisited sites onto a single curve. We present our results for 
p A (0) = 0.2 and use j ex = 1.45 and K = 0.86. 

The data collapse is excellent over both transient and 
asymptotic regimes, suggesting that the relation U(t) ~ 
K[p(t)\ 1 ' l:c holds to a high degree of accuracy for all times 
t. We attribute any deviations between the two curves 



7 



at late times to the onset of statistical fluctuations be- 
tween different runs of the data. Clearly we can improve 
the collapse at either early or late times by performing 
the linear regression in Fig. ^ on the appropriate part 
of the curve. We find that over the whole regime pre- 
sented in Fig. [h], r y ex ~ 1.45(1) works well. For brevity 
and clarity of presentation, we have presented our results 
for pa(0) = 0.2 only, although we have achieved equally 
good results for pa{0) = 0.3. 

The values j ex ~ 1.45 and 9 ex ~ 1.33 yield, using Eq. 
(0), 4>ab ^ 1-45 x 1.33 ~ 1.93. We therefore argue that, 
when infection is caused by both species, U(t) ~ t~^ AB 
with d>AB — 1-93. 



j ex ~ 1.45. An obvious goal for the future is to establish 
a theoretical framework for the exponents that occur in 
the exchange model and, in particular, for the observed 
simple relation between U(t) and P(t). Across the range 
of models we have studied, in both this and our previ- 
ous paper |57|, the fraction of uninfected walkers, U(t), 
and the fraction of unvisited sites, P(t), are seemingly 
related by U(t) ~ [P(t)] 7 . This remarkable relationship, 
which seems to hold reasonably accurately across the en- 
tire time regime accessible to simulation, merits further 
attention. In Table VI, we summarize the values of 7 for 
the various models we have studied (note, however, that 
there is some evidence [j57| that, for the first two models, 
the asymptotic value of 7 is smaller 1.39). 



IV. DISCUSSION AND SUMMARY 



In this paper, we have investigated the fraction of 
uninfected walkers U(t) in the ID q-state Potts model 
evolving at zero temperature. Our numerical results are 
consistent with previous work p6| , p4|j47| , suggesting that 
U(t) ~ with a nontrivial exponent <fi{q), where 

1/2 < <j)(q) < 3/2 for 2 < q < 00. Although the values of 
4>{q) for q — 2 and q — 00 are understood, the exact cal- 
culation of 4>{q) for all q, appears to be a hard problem. 
In analogy with the other models we have studied [p7| , we 
have reduced the familiar study of persistent sites Pit), 
to a limiting case in the study of uninfected walkers and 
noted that U{t) ~ [P(i)pW, with 4/3 < 7 < 3/2 for 
2 < q < 00. We discussed in Sec. I that, with regard to 
the persistence properties of the Potts model, it is pos- 
sible to probe the region 1 < q < 2. Although we have 
included a point for the state q = 1, where = = 0, 
in Fig. U and Fig. ||, we leave a study of the fraction of 
uninfected walkers in this latter regime as an open chal- 
lenge. 

We have also studied the ID exchange model, defined 
by the reactions A + A — > B, B + B — > A. Our numerical 
simulations indicate similarities between this diffusion- 
limited system and the Potts class of models, in that the 
particle density decays according to p(t) ~ t^ 1 ' 2 , inde- 
pendent of the initial particle density. The amplitude of 
the particle decay lies approximately halfway between the 
minimum and maximum values set by the Potts model 
suggesting that, within the context of the particle de- 
cay, the exchange model roughly mimics the behavior 
of two independent q = 4 Potts models. We have also 
shown that the fraction of sites unvisited by either an A 
or a B particle decays according to P(t) ~ t _e == where 
9 ex ~ 1.33. Our study of the fraction of uninfected walk- 
ers in the exchange model focused on two cases. In the 
first instance, we addressed the fraction of A walkers in- 
fected by only their own species and established that, in 
this case, U(t) ~ t~^ A where (j>A — 1.13. When cross 
species contamination was also included in our simula- 
tion, we found U(t) ~ t-^ AB with <\> A b ^ 1-93. We pre- 
sented numerical evidence that U(t) ~ K[p(t)] 1,l:L with 



Model 


7 


A%< — >QA 
A + B^V) 
g-state Potts Model 
exchange Model 


1.39 
1.39 

[4/3,3/2], q e [2, 00] 
1.45 



TABLE VI. Numerical results for the exponent 7 in various 
ID models for which U(t) ~ [P(i)] 7 (results for the first two 
models from ref. |57j). 

Finally, we note that all values of 7 obtained so far lie 
in the range 4/3 < q < 3/2 spanned by the g-state Potts 
model with 2 < q < 00. Up to now our analysis of 7 has 
been purely numerical, and naturally a more fundamen- 
tal understanding of this exponent is required. 
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